Deep Diversity: Extensive Variation in the Components of Complex Visual Systems across Animals.

Understanding the molecular underpinnings of the evolution of complex (multi-part) systems is a fundamental topic in biology. One unanswered question is to what the extent do similar or different genes and regulatory interactions underlie similar complex systems across species? Animal eyes and phototransduction (light detection) are outstanding systems to investigate this question because some of the genetics underlying these traits are well characterized in model organisms. However, comparative studies using non-model organisms are also necessary to understand the diversity and evolution of these traits. Here, we compare the characteristics of photoreceptor cells, opsins, and phototransduction cascades in diverse taxa, with a particular focus on cnidarians. In contrast to the common theme of deep homology, whereby similar traits develop mainly using homologous genes, comparisons of visual systems, especially in non-model organisms, are beginning to highlight a "deep diversity" of underlying components, illustrating how variation can underlie similar complex systems across taxa. Although using candidate genes from model organisms across diversity was a good starting point to understand the evolution of complex systems, unbiased genome-wide comparisons and subsequent functional validation will be necessary to uncover unique genes that comprise the complex systems of non-model groups to better understand biodiversity and its evolution.

Insect opsins and evo-devo: what have we learned in 25 years?

The visual pigments known as opsins are the primary molecular basis for colour vision in animals. Insects are among the most diverse of animal groups and their visual systems reflect a variety of life histories. The study of insect opsins in the fruit fly Drosophila melanogaster has led to major advances in the fields of neuroscience, development and evolution. In the last 25 years, research in D. melanogaster has improved our understanding of opsin genotype-phenotype relationships while comparative work in other insects has expanded our understanding of the evolution of insect eyes via gene duplication, coexpression and homologue switching. Even so, until recently, technology and sampling have limited our understanding of the fundamental mechanisms that evolution uses to shape the diversity of insect eyes. With the advent of genome editing and in vitro expression assays, the study of insect opsins is poised to reveal new frontiers in evolutionary biology, visual neuroscience, and animal behaviour. This article is part of the theme issue 'Understanding colour vision: molecular, physiological, neuronal and behavioural studies in arthropods'.

Multiple Mechanisms of Photoreceptor Spectral Tuning in Heliconius Butterflies.

The evolution of color vision is often studied through the lens of receptor gain relative to an ancestor with fewer spectral classes of photoreceptor. For instance, in Heliconius butterflies, a genus-specific UVRh opsin duplication led to the evolution of UV color discrimination in Heliconius erato females, a rare trait among butterflies. However, color vision evolution is not well understood in the context of loss. In Heliconius melpomene and Heliconius ismenius lineages, the UV2 receptor subtype has been lost, which limits female color vision in shorter wavelengths. Here, we compare the visual systems of butterflies that have either retained or lost the UV2 photoreceptor using intracellular recordings, ATAC-seq, and antibody staining. We identify several ways these butterflies modulate their color vision. In H. melpomene, chromatin reorganization has downregulated an otherwise intact UVRh2 gene, whereas in H. ismenius, pseudogenization has led to the truncation of UVRh2. In species that lack the UV2 receptor, the peak sensitivity of the remaining UV1 photoreceptor cell is shifted to longer wavelengths. Across Heliconius, we identify the widespread use of filtering pigments and co-expression of two opsins in the same photoreceptor cells. Multiple mechanisms of spectral tuning, including the molecular evolution of blue opsins, have led to the divergence of receptor sensitivities between species. The diversity of photoreceptor and ommatidial subtypes between species suggests that Heliconius visual systems are under varying selection pressures for color discrimination. Modulating the wavelengths of peak sensitivities of both the blue- and remaining UV-sensitive photoreceptor cells suggests that Heliconius species may have compensated for UV receptor loss.

Coordinated Gene Expression and Chromatin Regulation during Hydra Head Regeneration.

The cnidarian model organism Hydra has long been studied for its remarkable ability to regenerate its head, which is controlled by a head organizer located near the hypostome. The canonical Wnt pathway plays a central role in head organizer function during regeneration and during bud formation, which is the asexual mode of reproduction in Hydra. However, it is unclear how shared the developmental programs of head organizer genesis are in budding and regeneration. Time-series analysis of gene expression changes during head regeneration and budding revealed a set of 298 differentially expressed genes during the 48-h head regeneration and 72-h budding time courses. In order to understand the regulatory elements controlling Hydra head regeneration, we first identified 27,137 open-chromatin elements that are open in one or more sections of the organism body or regenerating tissue. We used histone modification ChIP-seq to identify 9,998 candidate proximal promoter and 3,018 candidate enhancer-like regions respectively. We show that a subset of these regulatory elements is dynamically remodeled during head regeneration and identify a set of transcription factor motifs that are enriched in the enhancer regions activated during head regeneration. Our results show that Hydra displays complex gene regulatory structures of developmentally dynamic enhancers, which suggests that the evolution of complex developmental enhancers predates the split of cnidarians and bilaterians.

A Bioinformatics Pipeline for Investigating Molecular Evolution and Gene Expression using RNA-seq.

Distilling and reporting large datasets, such as whole genome or transcriptome data, is often a daunting task. One way to break down results is to focus on one or more gene families that are significant to the organism and study. In this protocol, we outline bioinformatic steps to generate a phylogeny and to quantify the expression of genes of interest. Phylogenetic trees can give insight into how genes are evolving within and between species as well as reveal orthology. These results can be enhanced using RNA-seq data to compare the expression of these genes in different individuals or tissues. Studies of molecular evolution and expression can reveal modes of evolution and conservation of gene function between species. The characterization of a gene family can serve as a springboard for future studies and can highlight an important gene family in a new genome or transcriptome paper.

Molecular evolution and expression of opsin genes in Hydra vulgaris.

BACKGROUND: The evolution of opsin genes is of great interest because it can provide insight into the evolution of light detection and vision. An interesting group in which to study opsins is Cnidaria because it is a basal phylum sister to Bilateria with much visual diversity within the phylum. Hydra vulgaris (H. vulgaris) is a cnidarian with a plethora of genomic resources to characterize the opsin gene family. This eyeless cnidarian has a behavioral reaction to light, but it remains unknown which of its many opsins functions in light detection. Here, we used phylogenetics and RNA-seq to investigate the molecular evolution of opsin genes and their expression in H. vulgaris. We explored where opsin genes are located relative to each other in an improved genome assembly and where they belong in a cnidarian opsin phylogenetic tree. In addition, we used RNA-seq data from different tissues of the H. vulgaris adult body and different time points during regeneration and budding stages to gain insight into their potential functions. RESULTS: We identified 45 opsin genes in H. vulgaris, many of which were located near each other suggesting evolution by tandem duplications. Our phylogenetic tree of cnidarian opsin genes supported previous claims that they are evolving by lineage-specific duplications. We identified two H. vulgaris genes (HvOpA1 and HvOpB1) that fall outside of the two commonly determined Hydra groups; these genes possibly have a function in nematocytes and mucous gland cells respectively. We also found opsin genes that have similar expression patterns to phototransduction genes in H. vulgaris. We propose a H. vulgaris phototransduction cascade that has components of both ciliary and rhabdomeric cascades. CONCLUSIONS: This extensive study provides an in-depth look at the molecular evolution and expression of H. vulgaris opsin genes. The expression data that we have quantified can be used as a springboard for additional studies looking into the specific function of opsin genes in this species. Our phylogeny and expression data are valuable to investigations of opsin gene evolution and cnidarian biology.

Evolution of Phototransduction Genes in Lepidoptera.

Vision is underpinned by phototransduction, a signaling cascade that converts light energy into an electrical signal. Among insects, phototransduction is best understood in Drosophila melanogaster. Comparison of D. melanogaster against three insect species found several phototransduction gene gains and losses, however, lepidopterans were not examined. Diurnal butterflies and nocturnal moths occupy different light environments and have distinct eye morphologies, which might impact the expression of their phototransduction genes. Here we investigated: 1) how phototransduction genes vary in gene gain or loss between D. melanogaster and Lepidoptera, and 2) variations in phototransduction genes between moths and butterflies. To test our prediction of phototransduction differences due to distinct visual ecologies, we used insect reference genomes, phylogenetics, and moth and butterfly head RNA-Seq and transcriptome data. As expected, most phototransduction genes were conserved between D. melanogaster and Lepidoptera, with some exceptions. Notably, we found two lepidopteran opsins lacking a D. melanogaster ortholog. Using antibodies we found that one of these opsins, a candidate retinochrome, which we refer to as unclassified opsin (UnRh), is expressed in the crystalline cone cells and the pigment cells of the butterfly, Heliconius melpomene. Our results also show that butterflies express similar amounts of trp and trpl channel mRNAs, whereas moths express ∼50× less trp, a potential adaptation to darkness. Our findings suggest that while many single-copy D. melanogaster phototransduction genes are conserved in lepidopterans, phototransduction gene expression differences exist between moths and butterflies that may be linked to their visual light environment.

Hybrid Assembly of the Genome of the Entomopathogenic Nematode Steinernema carpocapsae Identifies the X-Chromosome.

Entomopathogenic nematodes from the genus Steinernema are lethal insect parasites that quickly kill their insect hosts with the help of their symbiotic bacteria. Steinernema carpocapsae is one of the most studied entomopathogens due to its broad lethality to diverse insect species and its effective commercial use as a biological control agent for insect pests, as well as a genetic model for studying parasitism, pathogenesis, and symbiosis. In this study, we used long-reads from the Pacific Biosciences platform and BioNano Genomics Irys system to assemble the most complete genome of the S. carpocapsae ALL strain to date, comprising 84.5 Mb in 16 scaffolds, with an N50 of 7.36 Mb. The largest scaffold, with 20.9 Mb, was identified as chromosome X based on sex-specific genome sequencing. The high level of contiguity allowed us to characterize gene density, repeat content, and GC content. RNA-seq data from 17 developmental stages, spanning from embryo to adult, were used to predict 30,957 gene models. Using this improved genome, we performed a macrosyntenic analysis to Caenorhabditis elegans and Pristionchus pacificus and found S. carpocapsae's chromosome X to be primarily orthologous to C. elegans' and P. pacificus' chromosome II and IV. We also investigated the expansion of protein families and gene expression differences between adult male and female stage nematodes. This new genome and more accurate set of annotations provide a foundation for additional comparative genomic and gene expression studies within the Steinernema clade and across the Nematoda phylum.

Item statistics derived from three-option versions of multiple-choice questions are usually as robust as four- or five-option versions: implications for exam design.

Different versions of multiple-choice exams were administered to an undergraduate class in human physiology as part of normal testing in the classroom. The goal was to evaluate whether the number of options (possible answers) per question influenced the effectiveness of this assessment. Three exams (each with three versions) were given to each of two sections during an academic quarter. All versions were equally long, with 30 questions: 10 questions with 3 options, 10 questions with 4, and 10 questions with 5 (always one correct answer plus distractors). Each question appeared in all three versions of an exam, with a different number of options in each version (three, four, or five). Discrimination (point biserial and upper-lower discrimination indexes) and difficulty were evaluated for each question. There was a small increase in difficulty (a lower average score on a question) when more options were provided. The upper-lower discrimination index indicated a small improvement in assessment of student learning with more options, although the point biserial did not. The total length of a question (number of words) was associated with a small increase in discrimination and difficulty, independent of the number of options. Quantitative questions were more likely to show an increase in discrimination with more options than nonquantitative questions, but this effect was very small. Therefore, for these testing conditions, there appears to be little advantage in providing more than three options per multiple-choice question, and there are disadvantages, such as needing more time for an exam.

Evolution of Sex-Biased Gene Expression and Dosage Compensation in the Eye and Brain of Heliconius Butterflies.

Differences in behavior and life history traits between females and males are the basis of divergent selective pressures between sexes. It has been suggested that a way for the two sexes to deal with different life history requirements is through sex-biased gene expression. In this study, we performed a comparative sex-biased gene expression analysis of the combined eye and brain transcriptome from five Heliconius species, H. charithonia, H. sara, H. erato, H. melpomene and H. doris, representing five of the main clades from the Heliconius phylogeny. We found that the degree of sexual dimorphism in gene expression is not conserved across Heliconius. Most of the sex-biased genes identified in each species are not sex-biased in any other, suggesting that sexual selection might have driven sexually dimorphic gene expression. Only three genes shared sex-biased expression across multiple species: ultraviolet opsin UVRh1 and orthologs of Drosophila Krüppel-homolog 1 and CG9492. We also observed that in some species female-biased genes have higher evolutionary rates, but in others, male-biased genes show the fastest rates when compared with unbiased genes, suggesting that selective forces driving sex-biased gene evolution in Heliconius act in a sex- and species-specific manner. Furthermore, we found dosage compensation in all the Heliconius tested, providing additional evidence for the conservation of dosage compensation across Lepidoptera. Finally, sex-biased genes are significantly enriched on the Z, a pattern that could be a result of sexually antagonistic selection.

Copy Number Variation and Expression Analysis Reveals a Nonorthologous Pinta Gene Family Member Involved in Butterfly Vision.

Vertebrate (cellular retinaldehyde-binding protein) and Drosophila (prolonged depolarization afterpotential is not apparent [PINTA]) proteins with a CRAL-TRIO domain transport retinal-based chromophores that bind to opsin proteins and are necessary for phototransduction. The CRAL-TRIO domain gene family is composed of genes that encode proteins with a common N-terminal structural domain. Although there is an expansion of this gene family in Lepidoptera, there is no lepidopteran ortholog of pinta. Further, the function of these genes in lepidopterans has not yet been established. Here, we explored the molecular evolution and expression of CRAL-TRIO domain genes in the butterfly Heliconius melpomene in order to identify a member of this gene family as a candidate chromophore transporter. We generated and searched a four tissue transcriptome and searched a reference genome for CRAL-TRIO domain genes. We expanded an insect CRAL-TRIO domain gene phylogeny to include H. melpomene and used 18 genomes from 4 subspecies to assess copy number variation. A transcriptome-wide differential expression analysis comparing four tissue types identified a CRAL-TRIO domain gene, Hme CTD31, upregulated in heads suggesting a potential role in vision for this CRAL-TRIO domain gene. RT-PCR and immunohistochemistry confirmed that Hme CTD31 and its protein product are expressed in the retina, specifically in primary and secondary pigment cells and in tracheal cells. Sequencing of eye protein extracts that fluoresce in the ultraviolet identified Hme CTD31 as a possible chromophore binding protein. Although we found several recent duplications and numerous copy number variants in CRAL-TRIO domain genes, we identified a single copy pinta paralog that likely binds the chromophore in butterflies.

Gene Duplication and Gene Expression Changes Play a Role in the Evolution of Candidate Pollen Feeding Genes in Heliconius Butterflies.

Heliconius possess a unique ability among butterflies to feed on pollen. Pollen feeding significantly extends their lifespan, and is thought to have been important to the diversification of the genus. We used RNA sequencing to examine feeding-related gene expression in the mouthparts of four species of Heliconius and one nonpollen feeding species, Eueides isabella We hypothesized that genes involved in morphology and protein metabolism might be upregulated in Heliconius because they have longer proboscides than Eueides, and because pollen contains more protein than nectar. Using de novo transcriptome assemblies, we tested these hypotheses by comparing gene expression in mouthparts against antennae and legs. We first looked for genes upregulated in mouthparts across all five species and discovered several hundred genes, many of which had functional annotations involving metabolism of proteins (cocoonase), lipids, and carbohydrates. We then looked specifically within Heliconius where we found eleven common upregulated genes with roles in morphology (CPR cuticle proteins), behavior (takeout-like), and metabolism (luciferase-like). Closer examination of these candidates revealed that cocoonase underwent several duplications along the lineage leading to heliconiine butterflies, including two Heliconius-specific duplications. Luciferase-like genes also underwent duplication within lepidopterans, and upregulation in Heliconius mouthparts. Reverse-transcription PCR confirmed that three cocoonases, a peptidase, and one luciferase-like gene are expressed in the proboscis with little to no expression in labial palps and salivary glands. Our results suggest pollen feeding, like other dietary specializations, was likely facilitated by adaptive expansions of preexisting genes-and that the butterfly proboscis is involved in digestive enzyme production.

Transcriptome-Wide Differential Gene Expression in Bicyclus anynana Butterflies: Female Vision-Related Genes Are More Plastic.

Vision is energetically costly to maintain. Consequently, over time many cave-adapted species downregulate the expression of vision genes or even lose their eyes and associated eye genes entirely. Alternatively, organisms that live in fluctuating environments, with different requirements for vision at different times, may evolve phenotypic plasticity for expression of vision genes. Here, we use a global transcriptomic and candidate gene approach to compare gene expression in the heads of a polyphenic butterfly. Bicyclus anynana have two seasonal forms that display sexual dimorphism and plasticity in eye morphology, and female-specific plasticity in opsin gene expression. Nonchoosy dry season females downregulate opsin expression, consistent with the high physiological cost of vision. To identify other genes associated with sexually dimorphic and seasonally plastic differences in vision, we analyzed RNA-sequencing data from whole head tissues. We identified two eye development genes (klarsicht and warts homologs) and an eye pigment biosynthesis gene (henna) differentially expressed between seasonal forms. By comparing sex-specific expression across seasonal forms, we found that klarsicht, warts, henna, and another eye development gene (domeless) were plastic in a female-specific manner. In a male-only analysis, white (w) was differentially expressed between seasonal forms. Reverse transcription polymerase chain reaction confirmed that warts and white are expressed in eyes only, whereas klarsicht, henna and domeless are expressed in both eyes and brain. We find that differential expression of eye development and eye pigment genes is associated with divergent eye phenotypes in B. anynana seasonal forms, and that there is a larger effect of season on female vision-related genes.

Genome Sequence of a Novel Iflavirus from mRNA Sequencing of the Butterfly Heliconius erato.

Here, we report the genome sequence of a novel iflavirus strain recovered from the neotropical butterfly Heliconius erato. The coding DNA sequence (CDS) of the iflavirus genome was 8,895 nucleotides in length, encoding a polyprotein that was 2,965 amino acids long.

Female behaviour drives expression and evolution of gustatory receptors in butterflies.

Secondary plant compounds are strong deterrents of insect oviposition and feeding, but may also be attractants for specialist herbivores. These insect-plant interactions are mediated by insect gustatory receptors (Grs) and olfactory receptors (Ors). An analysis of the reference genome of the butterfly Heliconius melpomene, which feeds on passion-flower vines (Passiflora spp.), together with whole-genome sequencing within the species and across the Heliconius phylogeny has permitted an unprecedented opportunity to study the patterns of gene duplication and copy-number variation (CNV) among these key sensory genes. We report in silico gene predictions of 73 Gr genes in the H. melpomene reference genome, including putative CO2, sugar, sugar alcohol, fructose, and bitter receptors. The majority of these Grs are the result of gene duplications since Heliconius shared a common ancestor with the monarch butterfly or the silkmoth. Among Grs but not Ors, CNVs are more common within species in those gene lineages that have also duplicated over this evolutionary time-scale, suggesting ongoing rapid gene family evolution. Deep sequencing (∼1 billion reads) of transcriptomes from proboscis and labial palps, antennae, and legs of adult H. melpomene males and females indicates that 67 of the predicted 73 Gr genes and 67 of the 70 predicted Or genes are expressed in these three tissues. Intriguingly, we find that one-third of all Grs show female-biased gene expression (n = 26) and nearly all of these (n = 21) are Heliconius-specific Grs. In fact, a significant excess of Grs that are expressed in female legs but not male legs are the result of recent gene duplication. This difference in Gr gene expression diversity between the sexes is accompanied by a striking sexual dimorphism in the abundance of gustatory sensilla on the forelegs of H. melpomene, suggesting that female oviposition behaviour drives the evolution of new gustatory receptors in butterfly genomes.